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We study the spin-half Heisenberg models on an anisotropic two-dimensional lattice which interpo- 
lates between the square-lattice at one end, a set of decoupled spin-chains on the other end, and the 
triangular-lattice Heisenberg model in between. By series expansions around two different dimer 
ground states and around various commensurate and incommensurate magnetically ordered states, 
we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly 
rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations and varying 
dimensionality. There is a large region of the usual 2-sublattice Neel phase, a 3-sublattice phase 
for the triangular-lattice model, a region of incommensurate magnetic order around the triangular- 
lattice model, and regions in parameter space where there is no magnetic order. We find that the 
incommensurate ordering wavevector is in general altered from its classical value by quantum fluc- 
tuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly 
critical. 



PACS Indices: 75.f0.-b., 75T0J., 75.40.Gb 

I. INTRODUCTION 

Quantum phases and phase transitions in two- 
dimensional Heisenberg models are of great current in- 
terest. JiJ Following much experimental, theoretical, and 
numerical work, a fairly comprehensive and consistent 
picture has emerged for the ordered phases of unfrus- 
trated Heisenberg models at zero and finite tempera- 
tures as well as for the quantum critical points sepa- 
rating them from the quantum disordered phases. [^]|| 
The case of frustrated Heisenberg models, potentially in- 
cludes much richer phenomena and is relatively less un- 
derstood. Here, the lattice models that have been stud- 
ied the most are the square-lattice nearest and second 
neighbor interaction (or J nn — J nnn ) model, where as a 
function of J nnn /J n „, the two and four sublattice colin- 
ear Neel phases are separated by a region of quantum 
disordered ground states. Q Theoretical arguments and 
various numerical results suggest that this intermediate 
phase is spontaneously dimerized. Other models that 
have received considerable attention are the triangular- 
lattice and Kagome-lattice Heisenberg models. ||-|7J The 
former appears to be magnetically ordered in the three- 
sublattice non-colinear pattern, whereas the latter almost 
surely has a non-magnetic ground state. The nature of 
the magnetically disordered phases and low energy ex- 
citations in these frustrated Heisenberg models are sub- 
jects of considerable interest. 

Here we study a class of Heisenberg models that inter- 



polate between the square-lattice, the triangular-lattice 
and the linear chain Heisenberg models. It can be de- 
fined in terms of a square-lattice of spin-half operators, 
with a nearest neighbor Heisenberg interaction Ji and a 
second neighbor interaction J\ along only one of the di- 
agonals of the lattice. In the limit J\ — > 0, we recover 
the square-lattice Heisenberg model. In the limit J2 — > 0, 
the model reduces to decoupled spin-chains running along 
the diagonals of the square-lattice. For small J2/J1 the 
model consists of weakly coupled chains with frustrated 
interchain coupling. For J2 = J%, the model is exactly 
equivalent to the triangular-lattice Heisenberg model. 

This model is of particular interest for at least two 
reasons. First, within a single Hamiltonian, it provides 
a way of interpolating between several well-known one 
and two dimensional models. One can study the role 
of frustration in going from one to two dimensions, (^] 
as well as the stability of commensurate and incommen- 
surate magnetic order with respect to quantum fluctu- 
ations. [|| Second, the model is of direct relevance to 
the magnetic phases of various quasi-two-dimensional or- 
ganic superconductors. It has been argued that this 
model should describe the spin degrees of freedom of 
the insulating phase of the layered molecular crystals, 
k-(BEDT-TTF) 2 X. (To) These materials probably have 
J1/J2 ~ 0.3 — 1. This ratio varies with the anion X 
and should vary with uniaxial stress applied within the 
layer along the diagonal. JllJ] A very slight variant on this 
model with J1/J2 ~ 4 has been proposed to describe the 
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insulating phase of (9-(BEDT-TTF) 2 RbZn(SCN) 4 which 
has a spin gap. Recent studies of the associated Hubbard 
model on an anisotropic triangular lattice at half-filling 
suggest that the wave vector associated with the spin ex- 
citations determines the type of superconducting order. 
p2[ Consequently, it appears that a good understand- 
ing of this Heisenberg model may be a pre-requisite to 
understanding organic superconductivity. 

In this paper we study the zero-temperature phase di- 
agram of this model by series expansions around various 
dimerized and magnetically ordered phases. Because of 
frustration, and large coordination number of the lattice, 
it would be difficult to get a comparable treatment of this 
model by any other numerical method. We will only con- 
sider Ji > 0, Ji > 0, that is both couplings are antifer- 
romagnetic, which introduces frustration into the model. 
We find that the Neel order of the square-lattice Heisen- 
berg model persists up to J1/J2 ^5 0.7, where it goes 
continuously to zero. It is interesting to note that Neel 
order exists beyond it regime of classical stability, provid- 
ing an example of quantum order by disorder. |9) In the 
region 0.7 < J1/J2 ^ 0.9, there is no magnetic order and 
the ground-state is dimerized in the columnar pattern. 
For larger values of J1/J2, there is incommensurate or 
spiral long-range order. Except for the triangular-lattice 
point (J1/J2 = 1), we find that the ordering wavevector 
is different from the classical Heisenberg model. At the 
triangular-lattice point, the dimer expansions indicate a 
magnetic instability at an ordering wavevector which co- 
incides with the classical value. Furthermore, our Ising 
expansion series in this case reduce to previously derived 
series expansions directly for the triangular-lattice model. 

§ 

The region of large J1/J2 is very interesting and is 
not fully resolved by our series expansions. In the limit 
of small J2, we have Heisenberg spin-chains coupled by 
frustrating zig-zag interactions. For just two such chains, 
this problem was studied by Affleck and White. |l3| This 
zig-zag ladder has incommensurate spin correlations, a 
small gap and a spontaneous dimerization, which breaks 
the degeneracy of the weak zig-zag bonds but leaves the 
strong chain bonds equivalent. Our numerical results 
suggest that such a state is definitely not favored for 
the two-dimensional system. This is not hard to under- 
stand as such a broken symmetry phase can only have 
the strong bonds between every alternate pair of chains. 
We suggest three possible scenarios for this phase: (i) 
It has weak incommensurate spiral order, which vanishes 
rapidly (with a power greater than unity) as J2 — ► 0, 
(ii) This phase is dimerized along the strong J\ bonds 
and has a spin-gap, (iii) The phase remains critical with- 
out any long-range order. The ground state energy ob- 
tained from the Ising expansions around a spiral phase 
and around a dimer state along the strong bonds remain 
nearly equal through much of this phase. In all three 
cases the system is close to being critical. 



The elementary excitations of this phase and how they 
are connected to the excitations in isolated spin-chains 
should be very interesting to investigate. Q The low- 
energy properties of models with non-collinear order are 
described by the SO (3) non-linear sigma model. |L5 16 
In the quantum disordered phase of this model the ele- 
mentary excitations are deconfined spin-1/2 objects, or 
spinons, in contrast to the magnons or spin-one 

triplets that are the elementary excitations in most two- 
dimensional antiferromagnets. Thus, this model, with 
spiral classical order provides us with candidate systems 
to look for spinon excitations. 

The plan of the paper is as follows. In Sec. II we 
discuss the model and the various series expansions that 
are carried out. In Sec. Ill we discuss the results for the 
phase diagram, the ground state energy, the excitation 
spectra and the correlation functions of the model. Fi- 
nally, in the last section we present our conclusions and 
discuss some future directions. 
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FIG. 1. The anisotropic triangular lattice with coupling Ji 
and J2. 



II. SERIES EXPANSIONS 

The Heisenberg antiferromagnet on an anisotropic tri- 
angular lattice is equivalent to the Heisenberg antiferro- 
magnet on a square lattice with nearest neighbor inter- 
actions and one of diagonal next-nearest neighbor inter- 
actions as shown in Fig. 1. The Hamiltonian is: 



H 



Ji S,j • S n + J 2 Si • Sj 



(1) 



(in) 



where {ij} means pairs of nearest neighbors and (in) 
means one of diagonal next-nearest neighbor pairs, as 
illustrated in Fig. 1. We denote the ratio of the cou- 
plings as y, that is, y = J\/ 'J2. In the present paper, we 
study only the case of antifcrromagnetic coupling, where 
both Ji and J 2 are positive. In the limit Ji = 0, the 
model is equivalent to decoupled one-dimensional Heisen- 
berg chains. In the limit J\ = 0, the model is equivalent 
to the two-dimensional square-lattice Heisenberg model. 
For Ji — J 2 , the model is equivalent to the Heisenberg 
model on a triangular lattice. 

Let's discuss the classical ground state of this system: 
If we assume (Si) lies in the x — z plane, we can rotate 
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all the spins related to a reference spin so we have a 
ferromagnetic ground state. The classical ground state 
of this system has nearest neighbor spins differ by an 
angle q, and spin joined by a diagonal bond differ by 2q, 
where q is 



q = 



7T, Ji < J 2 /2; 

arccos(— J2/2J1), J\ > J2/2. 



(2) 



The angle q as a function of Ji/(Ji + J2) is shown as the 
solid curve in Fig. g. 

After this rotation, the transformed Hamiltonian is 



H — Hi + J1H2 + J2H3, 



(3) 



where 



A. Ising Expansions 

To construct a T = expansion about the classical 
ground state for this system, one has to introduce an 
anisotropy parameter A, and write the Hamiltonian for 
the Heisenberg-Ising model as 



H = H +xV , 



(5) 



where 



Hq — Hi - t S* , 

i 

V= JiH 2 + J 2 H 3 + tY,S* 



(6) 



The last term of strength t in both Hq and V is a lo- 
cal field term, which can be included to improve con- 
vergence. The limits x = and x = 1 correspond to 
the Ising model and the isotropic Heisenberg model, re- 
spectively. The operator Hq is taken as the unperturbed 
Hamiltonian, with the unperturbed ground state being 
the usual ferromagnetically ordered state. The operator 
V is treated as a perturbation. 



Hi = Ji cos(2g) S?S* + J 2 cos(g) ^ S!f , 

{in) (ij) 

H 2 = J2 S!Sl + cos(2g)Sf S* n + sin(2 g )(5f S% - S*S z n ) (4) 
(in) 

H 3 = S? + co S (q)SfS* + smbKS'Sf ~ SfS?) . 

We have studied this system by using linked-cluster 
expansion methods at zero-temperature including Ising 
expansions, and dimer expansions about two different 
dimerization patterns. The linked-cluster expansion 
method has been previously reviewed in several articles 
pTsj po[ , and will not be repeated here. Here we will only 
summarize the expansion methods used, and the results 
derived from them will be presented. The series coeffi- 
cients are not given here, but are available upon request. 
Previous studies using these methods have agreed quan- 
titatively with other numerical studies (including quan- 
tum Monte Carlo and exact diagonalization) for unfrus- 
trated spin models. For frustrated two-dimensional mod- 
els, both the Monte Carlo and the exact diagonalization 
methods face problems (due to the minus signs in case 
of the former and tremendous variations with size and 
shapes of the clusters in case of the latter). Thus, we 
believe, it would be difficult to get a comparable numer- 
ical treatment of this strongly frustrated model by other 
numerical methods. 
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FIG. 2. Phase diagram of the model. The wavevector q 
is plotted as a function of Ji/(Ji + J2). The solid curve is 
the classical ordering wavevector, while the points are esti- 
mates of the location of the minimum triplet gap from both 
columnar dimer expansions (solid points) and diagonal dimer 
expansions (open points). The long dashed lines are the fit 
given in Eq. (y_6(). The vertical short dashed lines indicate the 
region of spontaneously dimerized phase, and vertical dotted 
line indicate another possible transition point. 

To determine whether the system is in the ordered 
phase, we need to compute the order parameter M de- 
fined by 



M = ^E(S!) 



(7) 



where the brackets denote ground-state expectation val- 
ues. Actually, in this expansion, we can choose q different 
from that given in Eq. (0), the favored q should give a 
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lower energy. If q = tt, we are expanding about the usual 
Neel phase. 

Ising series have been calculated for the ground state 
energy per site, Eq/N, the order parameter M, and the 
triplet excitation spectrum (for the case of q = tt only) 
for several ratios of couplings x and (simultaneously) for 
several values of t up to order a; 10 . This calculation in- 
volves 116137 linked-clusters of up to 10 sites. The series 
for the case of y = 1 have been computed previously || , 
and our results agree with these previous results. This is 
a highly non-trivial check on the correctness of the expan- 
sion coefficients, as the set of (anisotropic square-lattice) 
graphs considered here are completely different from the 
(triangular-lattice) graphs considered in Ref. || and the 
representation of the Hamiltonian as a Heisenberg-Ising 
model is also superficially different. 

At the next stage of the analysis, we try to extrapolate 
the series to the isotropic point [x = 1) for those values 
of the exchange coupling parameters which lie within the 
ordered phase at x — 1. For this purpose, we first trans- 
form the series to a new variable 



l-(l-ar) 



1/2 



(8) 



to remove the singularity at x — 1 predicted by the spin- 
wave theory. This was first proposed by Huse |2lJ and 
was also used in earlier work on the square lattice case 
p2| . We then use both integrated first-order inhomo- 
geneous differential approximants p3fl and Pade approx- 
imants to extrapolate the series to the isotropic point 
8 = 1 (x = 1). The results of the Ising expansions will 
be presented in later sections. 



and the other is the diagonal covering displayed in Fig. 
|(a). 



1. Columnar dimer expansions 

In this expansion, we take the columnar bonds (D) 
denoted by the thick solid bonds in Fig. || as the un- 
perturbed Hamiltonian Hq, and the rest of the bonds as 
perturbation, that is, we define the following "dimerized 
Ji — J2 model" 

H = J 2 S, ; • Sj + x J 2 22 Si ■ Sj + J\ 22 ^ ' Sj 

(9) 

We take the first term in above H as the unperturbed 
Hamiltonian, and the second term as perturbation. The 
unperturbed ground state is a product state of near- 
est neighbor singlet dimers and the perturbation couples 
these among themselves and with the pair-triplet states. 
At x = 1, one recovers the original Hamiltonian (|l|). 

Two types of series can be obtained in this way. By 
fixing y, we can compute series in the single variable x 
up to order L, 



(10) 



i=0 



In the other approach we keep both x and y as expansion 
parameters and obtain double series of the form 
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FIG. 3. Columnar dimerization pattern, the thick solid 
lines represent strongest bonds J2, and thin solid lines and 
dashed lines represent weaker bonds XJ2 and xJi, respec- 
tively. Note that the dimer centers lie on a rectangle lattice 
with spacing 2a in x-direction and a in y-direction. In char- 
acterizing the excitation spectrum we set the lattice spacing 
in both the x- and y-direction to be 1. 



B. Dimer Expansions 

For this system, we have carried out two types of dimer 
expansions corresponding to two types of dimer cover- 
ings. One is the columnar covering displayed in Fig. 0, 



= ££*. 

i=0 j=0 



(ii) 



where the coefficients dij are computed up to order i = L. 

The series has been computed up to order L = 9 for the 
ground-state energy Eq, and up to order L = 8 for the 
triplet excitation spectrum A(k Xl k y ). This calculation 
involves 49684 linked clusters of up to 9 sites. 

Note that if the centers of the dimers are treated as 
lattice sites, one obtains a rectangular lattice with spac- 
ing 2a in ^-direction and a in y-direction, where a is the 
lattice spacing of original lattice. In characterizing the 
excitation spectrum presented in the later sections we 
set the lattice spacing in both x- and y-directions to be 
1 . In this notation, the wavevector corresponding to the 
classical ordering is (q ~ tt, q), which varies from (0, tt) to 
(— 7r/2,7r/2) as J1/J2 increases from to 00. 

The leading order series for the triplet excitation spec- 
trum is: 

A(k x ,k y )/J 2 = 1 + x[~\ cos(k x ) + (1 - Ji/2 J 2 ) cos(A; y ) 



-(Ji/2J 2 ) cos(k x + k y )} + 0(x 2 ) 



(12) 
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Note that the minimum of the above spectrum may not 
coincide with the ordering wavevector (g — tt, q) of the 
classical system. 



xJ, 



(a) 




FIG. 4. (a) Diagonal dimerization pattern, the thick 
dashed lines represent strongest bonds Ji, and thin dashed 
lines and solid lines represent weaker bonds xJi and £ J2, re- 
spectively, (b) The equivalent lattice to (a) used in series 
expansions. Again note that the dimer centers lie on a rectan- 
gle lattice with spacing 2a in ^-direction and a in y-direction. 
In characterizing the excitation spectrum we set the lattice 
spacing in both x- and {/-direction to be 1. 



2. Diagonal dimer expansions 

In this expansion, we take the diagonal bonds (D) de- 
noted by the thick dashed bonds in Fig. |](a) as unper- 
turbed Hamiltonian Hq, and rest of the bonds as pertur- 
bation. That is wc define the following "dimerized J\ — J2 
model" 

H = J\ Si-Sj+x J\ ^2 • Sj + J% Si ■ Sj 

{i,ri)eD [i,n)$.D (ij) 

(13) 

Again at x = 1, one recovers the original Hamiltonian 

For convenience, we have transformed the lattice 
shown in Fig. ^(a) to the equivalent one shown in Fig. 
^(b). Again the dimer centers in Fig. [|(b) lie on a rect- 
angular lattice with spacing 2a in x-direction and a in y- 
direction. In characterizing the excitation spectrum pre- 
sented in the later sections we also set the lattice spacing 
in both x- and y-directions to be 1. In this representation 
the wavevector corresponding to the classical ordering is 
(tt — 2q,q), which varies from (— ir,n) to (0,7r/2) as J\j J% 
increases from to 00. 



As in the case of columnar dimer expansions, two types 
of series have been computed up to order x for the 
ground-state energy E , and up to order a; 8 for the triplet 
excitation spectrum A(k x ,k y ). This calculation involves 
the same linked clusters as the columnar dimer expan- 
sions. 

The leading order series for triplet excitation spectrum 

is 

A(k x ,k y )= 1 + cos(k x ) + ( J 2 /2 Ji) cos(fc y ) 

-(J 2 /2Ji) cos{k x + k y )} + 0{x 2 ) (14) 

Note that the minimum of A is at (ir — 2q 1 q), with q 
defined in Eq. (^|), that is just the ordering wavevector 
of the classical system. 



III. RESULTS 

Having obtained the series for the various expansions 
above we present in this section the results of series anal- 
ysis. We use integrated first-order inhomogeneous dif- 
ferential approximants |p3] | and Pade approximants to 
extrapolate the series. 



-0.25 



-0.3 



z 

\ 

o 
H 



-0.35 - 



-0.4 - 




FIG. 5. The ground-state energy per site Eo/N(Ji +J2) as 
function of Ji/(Ji + J2), obtained from the Ising expansions 
about the Neel order (i.e. q = tt) (A), the Ising expansions 
about the classical spiral order (i.e. q 7^ tt defined by Eq. 
(^) (B), the Ising expansions using the q obtained by dimer 
expansions (C), the columnar dimer expansions (D) and the 
diagonal dimer expansion (E). The result at J2 = is the ex- 
act results for the ground state energy of a antiferromagnetic 
spin-chain [28]. 
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A. Ground state energy 

Fig. H shows the ground state energy per site E /N 
for the Hamiltonian in Eq. ([!]) obtained from various 
expansions, We see that the Ising expansions about the 
Neel order give the lowest energy and also show the best 
convergence in the region of J\j Ji < 0.7. The Ising ex- 
pansions about the classical ground state (i.e. q tt 
defined by Eq. (^)) gives the lowest energy in the re- 
gion of 0.9 < J1/J2 ^ 1.2, the columnar dimer expan- 
sions gives the lowest energy for 0.65 < J1/J2 0.9, and 
the diagonal dimer expansions give the lowest energy for 
J1/J2 ^ 1-2. Whenever the dimer expansions indicate 
an ordering wavevector different from the classical one, 
we have also computed the series expansion around the 
spiral-state with this modified wavevector. These always 
give lower ground state energies than those for the clas- 
sical wavevector. Also, for J1/J2 ^ 1.2, these modified 
Ising expansions lead to ground state energies that are 
almost identical with the diagonal dimer expansions. 
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theory. |24| The staggered magnetization associated with 
the Neel order exists up to J\j Ji — 0.7, rather than up 
to J1/J2 = 0.5, as in the classical system. We can also 
see that the magnetisation M vanishes over the region 
0.7 < J1/J2 ;$ 0.9, suggesting that this may belong to a 
quantum disordered phase. For J1/J2 ^ 0.9, the system 
may have spiral order. 

The fact that the Neel order survives beyond the classi- 
caly stable region, is an example of promotion of collincar 
order by quantum fluctuations. A similar stabilization 
was found in the 1/4 depleted square-lattice Heisenberg 
model relevant to the material CaVaOy by Kontani et 
al. [^H This can be viewed as a part of the more gen- 
eral phenomena of quantum order by disorder, || where 
quantum fluctuations help select and stabilize an appro- 
priate type of order, typically collinear, in face of near 
classical degeneracy. 




FIG. 7. Plot of triplet excitation spectra A(k x ,k y ) (de- 
rived from the Ising expansions about Neel order) along 
high-symmetry cuts through the Brillouin zone for coupling 
ratios J1/J2 = 0,0.3,0.6. 



FIG. 6. The order parameter M versus J1/J2. The solid 
points with error bars are the estimates from the Ising expan- 
sions about Neel order (i.e. q = tt) , while the open points are 
the estimates from Ising expansion about the classical spiral 
order (i.e., q defined by Eq. (0)). 



B. Magnetisation 

The order parameter M obtained from the Ising ex- 
pansions about the classical ground state is shown in Fig. 
||. The dependence on the amount of magnetic frustra- 
tion is very similar to that obtained from linear spin wave 
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FIG. 8. Plot of triplet excitation spectrum A(k x ,k y ) for 
the columnar dimerized Ji — J2 model at x — 0.5 for cou- 
pling ratios J1/J2 = 0.6,0.8,0.9,1,1.1,1.25. The big circles 
indicate the locations of classical ordering wavevectors. 

C. Triplet excitation spectrum 

Fig. ^ shows the triplet excitation spectra A(k x ,k y ) 
for J1/J2 = 0,0.3,0.6 obtained from the Ising expan- 
sions about Need order, where we can see the energy gap 
vanishes at (71-, 71") and its equivalent point. 

We now present the spectra calculated from the colum- 
nar dimer expansions. We have previously noted that 
the leading order columnar dimer expansions give the 
minimum of the triplet excitation spectrum close to the 
classical ordering wavevector. Fig. |8| shows the triplet 
excitation spectra A(k x ,k y ) for the columnar dimerized 
Ji — J2 model at x = 0.5 and coupling ratios J1/J2 = 
0.6, 0.8, 0.9, 1, 1.1, 1.25. It can be seen that the minimum 
gap is located at (0, n) for J1/J2 ^ 0.7, or (0,0) for 
J1/J2 Z 1-25. For J1/J2 over the region 0.8 - 1.1, the 
minimum is close to the classical ordering wavevector. 



FIG. 9. Plot of triplet excitation spectrum A(k x ,k y ) de- 
rived from the columnar dimer expansions for the coupling 
ratios J1/J2 = 0.75, 0.85 and x = 1. The big full circles indi- 
cate the locations of classical ordering wavevector. 

Extrapolating the columnar dimer series to x — 1 by 
the integrated differential approximants, one can get the 
triplet excitation spectra for the fairly narrow parameter 
range 0.7 < J1/J2 ^ 0.9, and the results for Ji/ J 2 = 0.75 
and 0.85, are shown in Fig. ^. We find that the mini- 
mum gap is nonzero for these couplings. The minimum 
triplet gap as a function of J1/J2 are given in Fig. [lC| . 
We find that up to J1/J2 — 0.75, the minimum gap is 
located at (0,7r), and for J1/J2 ^ 0.75, the system is in 
an incommensurate phase, with the location of minimum 
gap given in Fig. ^|. 

The excitation spectra for the regime of small J2/J1 
will be discussed in a later section. 
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FIG. 10. The minimum triplet excitation gap A vs J1/J2 
derived Irom the columnar dimer expansions. The full points 
indicate the minimum gap is located at (0, n). 



D. Phase Transitions and Critical Exponents 

In this section we describe various estimations of the 
phase boundaries and the calculations of critical expo- 
nents associated with the transitions. We note that the 
estimates of critical exponents are not particularly ac- 
curate. The uncertainties shown are a measure of the 
spread among the Pade approximants. Given the short 
length of the series, they are best regarded as effective 
exponents. 

First, combining the results of the Ising and the colum- 
nar dimer expansions, we estimate that Neel order dis- 
appears at J1/J2 = 0.68(3). The phase adjacent to that 
appears to be spontaneously dimerized in the columnar 
pattern. This phase is analogous to the middle spon- 
taneously dimerized phase in J nn — J nnn square lattice 
model. Q| Our results are consistent with the transition 
from the Neel to the dimerized phase being a second- 
order transition, since, near the critical point, the ground 
state energy from the columnar dimer expansion matches 
very smoothly those from the Ising expansions. For a 
first-order transition, there should be a discontinuity in 
the slope at the transition. Similarly, we estimate the 
transition point between spontaneously dimerized phase 
and the spiral phase to be at J\j J% = 0.91(3). 

Phase boundaries for the columnar dimerized Hamil- 
tonian can be more fully studied by the columnar dimer 
expansions. By applying the Dlog Pade approximants 
to the triplet series at the wavevector of the minimum 



gap, we can locate the critical points at which the gap 
closes. The results are presented in Fig. O, where in 
our analysis we have assumed that the minimum gap is 
located at the classical ordering wavevector for J\ j J2 in 
the range 0.8 — 1.2. The associated critical exponent v 
is about 0.8(1), with not much variation with J\/J2- It 
would be interesting to calculate the spectral weight for 
the triplet excitations to see if they vanish, suggesting 
the appearance of spinon excitations. |l4j These results 
are not accurate enough to decide if the transitions be- 
long to the 0(3) nonlinear sigma model universality class 
(y = 0.70) (2(| or the 0(4) exponents associated with 
the SO(3) nonlinear sigma model (v — 0.74) |]l5| , ^7| . For 
J1/J2 in the region 0.7 — 0.9, we find critical points x c 
larger than unity. This is evidence that this phase is 
spontaneously dimerized. 




0.4 0.6 
J 1 /(J 1 +J 2 ) 



FIG. 11. The critical point x c as a function of Ji/(Ji + J2) 
as obtained from Dlog Pade approximants to the minimum 
energy gap of columnar dimer expansions. The locations of 
minimum gap are (0,7r) for J1/J2 < 0.6, or near classical 
ordering wavevector (q — n,q) for 0.95 < J1/J2 ^ 1-2, or (0,0) 
for J1/J2 > 1.25. 

For diagonal dimer expansions, we can also get the 
phase boundary by using the Dlog Pade approximants to 
locate the critical points where the gap vanishes. The 
results are shown in Fig. 12. The associated critical 



exponent v is found to be approximately 0.75(10), and 
also does not vary much with J^l J\- For J2/J1 1-75 
the minimum of the energy gap is clearly (— tt, 7r), which 
is the classical ordering wavevector. 
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E. Weakly Coupled Zig-Zag Chains 



(0,0) 



(f-2q,q) 



/ 



/ 

/ (-tt.tt) 



Diagonal dimer expansions 



0.4 0.6 

V(J,+J a ) 



FIG. 12. The critical point function of Ji/(Ji + J2) 

as obtained from Dlog Pade approximants to the minimum 
energy gap of diagonal dimer expansions. The locations of 
minimum gap are (0,0) for J2/J1 < 0.75, or near classical 
ordering wavevector (ir — 2q,q) for 0.75 < J2/J1 J$ 1.75, or 
(-7r,7r) for J2/J1 > 1.75. 
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FIG. 13. Plot of triplet excitation spectrum A(fc a; ,/cj / ) 
along selected contours in Brillouin zone (derived from the di- 
agonal dimer expansions) for coupling ratios J2/J1 = 0. The 
big full circle indicates the position of the classical ordering 
wavevector and q x = n(n — 2q)/q = for J2/J1 = 0. 



In this section we discuss the excitation spectra for 
small J2/J1, when the system can be thought of as weakly 
coupled spin-chains. The coupling between the chains is 
frustrated, and has been called zig-zag coupling in the 
literature. 



J 2 /J,=0.2 

■r 

/ 



(0.0) (O.ti/2) (-".") (0.0) 

(k,, k ( ) 



FIG. 14. Plot of triplet excitation spectrum A(k x ,k y ) 
along selected contours in Brillouin zone (derived from the 
diagonal dimer expansions) for coupling ratios J2/J1 = 0.2. 
The big full circle indicates the position of the classical or- 
dering wavevector and q x — n(n — 2q)/q = — 0.127T for 
J2/J1 = 0.2. 

For small values of J2/J1, the critical points in the 
diagonal dimer expansions lie very close to unity, and 
it is very difficult to locate the wavevector of the gap 
minimum at x = 1. Fig. p"3| gives the dispersion at x — 1 
for the case J2 = (i.e., decoupled spin chains). Here we 
get a small but nonzero minimum gap, due to the fact 
that we have not taken into account the singular behavior 
of the spin-chain. This can be done by extrapolating the 
series in a new variable 



5 = 1 - (1 - xf/\ 



(15) 



in which case the gap does vanish. For J2 non-zero, we 
do not know if x = 1 is gapped or critical, and hence 
do not know which extrapolation variable to use. For 
simplicity, we discuss primarily an extrapolation in the 
original variable, which assumes there is no singularity at 
x = 1. Note, however, that even with this extrapolation 
we do get a flat dispersion between (0,0) and (0, 7r/2), 
as expected. As we increase J2/J1 (up to J2/J1 — 0.2), 
we still get a very flat dispersion over the region (0, 0) 
and (0, 7r/2). The gap over this region is also very small 
(similar to that for J2 =0), and in fact it is zero within 
errorbars if we extrapolate the series using the variable 
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S. As an example, the dispersion is shown in Fig. |l4| for 
J2/J1 — 0.2 (extrapolated assuming no critical points 
for x < 1). At larger values of J2/J1, for example for 
J2/J1 = 0.5, the dispersion is more pronounced and is 
shown in Fig. [l5|. We can see that the position of the 
minimum gap is not located exactly at the classical or- 
dering wavevector, but deviates from it along the direc- 
tion connecting (0, 7r/2) and (— 7r,7r). This location of 
the minimum gap is plotted in Fig. 0. 




1 — ' — ' — ' — ' — 1 — ' — ' — ' — ' — 1 — ' — ' — ' — ' — 1 — 1 — 1 — ' — 1 — 1 

(0,0) (O.Ti/Z) (-it.t) (0,0) (q„,7T) 

(K k,) 



FIG. 15. Plot of triplet excitation spectrum A(k x ,k y ) 
along selected contours in Brillouin zone (derived from the 
diagonal dimer expansions) for coupling ratios J2/J1 = 0.5. 
The big full point indicates the position of classical ordering 
wavevector, while the big open square indicates the position 
of minimum gap. The solid (open) points in rightmost win- 
dows are spectra along contour connecting (0, 0) and the clas- 
sical ordering wavevector (the position of minimum gap) by 
a straight line up to k y = n. 

The wavevector q obtained from both columnar dimer 
expansions and diagonal dimer expansions can be fitted 
very well by the simple form: 

_ f arccos[-(3J 2 /Ji - 2)/2], 0.75 > J x /J 2 < 1; 
q ~ { arccos{-[2(3Ji/J 2 - 2)]" 1 }, J x /J 2 > 1. 

(16) 

which is shown by a dashed lined in Fig. ^. 

Using the above results for q, we can obtain more reli- 
able results for the ground state energy by the Ising ex- 
pansion, and the results are shown in Fig. ^| The results 
for the magnetisation M obtained using this wavevector 
agree (within errorbars) with the results obtained using 
the classical wavevector, shown in Fig. |^ and, hence, are 
not shown. 



IV. CONCLUSIONS 

In this paper we have studied a class of Heisenberg 
models, which interpolate between the square-lattice, the 
triangular-lattice and the linear-chain models. As dis- 
cussed in the introduction, these models should describe 
the magnetic properties of the insulating phase of a range 
of layered organic superconductors. Furthermore, using 
uniaxial stress it may be possible to vary the parameter 
J\ I Ji and induce some of the quantum phase transitions 
considered in this paper. We find that the model has a 
particulary rich phase diagram: Neel, spiral, dimerized, 
and possible critical phases occur. Except, perhaps, for 
the region of weakly coupled chains the phase diagram 
is rather well determined by our series expansions, and 
it would be difficult for any other numerical method to 
have comparable accuracy. Some of the most notable 
features are that the two-sublattice Neel phase appears 
stable beyond its classical regime of stability and in the 
spiral phase, the ordering wavevector is different from 
the classical ordering wavevector except for the case of 
the triangular-lattice Heisenberg model. The region of 
weakly coupled chains is nearly critical. In studying the 
instability of the dimerized phase to the spiral phase, the 
critical exponents for various dimer expansions give val- 
ues of v in the range 0.65 < v < 0.9 which is consistent 
with either 0(4) or 0(3) universality class. This, as well 
as the possibility of spinon excitations in the quantum 
disordered phases of this model, deserve further atten- 
tion. 

ACKNOWLEDGMENTS 

This work is supported in part by a grant from the 
National Science Foundation (DMR-9616574) (R.R.P.S), 
the Gordon Godfrey Bequest for Theoretical Physics at 
the University of NSW, and by the Australian Research 
Council (ZW. and R.M.). The computation has been 
performed on Silicon Graphics Power Challenge and Con- 
vex machines. We thank the New South Wales Centre 
for Parallel Computing for facilities and assistance with 
the calculations. We thank J. Merino, J. Oitmaa, C.J. 
Hamer and S. Sachdev for helpful discussions. We thank 
C.J. Hamer and J. Oitmaa for providing some of the com- 
puter codes we used. 



e-mail address: w.zheng@unsw.edu.au 
e-mail address: ross@newt.unsw.edu.au 
e-mail address: rrpsing@ucdavis.edu 



10 



[1] 

[2] 
[3] 

[4] 

[5] 

[6] 
[7] 
[8] 

[9] 

[10] 

[11] 
[12] 

[13] 
[14] 



S. Chakravarty, B. I. Halperin and J. Ye, Phys. Rev. B 
39, 2344 (1989); A. V. Chubukov, S. Sachdev and J. Ye, 
Phys. Rev. B 49, 11919 (1994). 

A. Auerbach, Interacting Electrons and Quantum Mag- 
netism, (Springer - Verlag, New York, 1994). 
S. Sachdev, in Low Dimensional Quantum Field Theo- 
ries for Condensed Matter Physicists, edited by Y. Lu, 
S. Lundqvist, and G. Morandi, (World Scientific, Singa- 
pore, 1995). 

H. J. Schulz, T. A. L. Ziman, and D. Poilblanc, J. Phys. 
I France 6, 675 (1996); J. Oitmaa and Z. Weihong, 
Phys. Rev. B 54, 3022 (1996), and references therein; O. 
Sushkov, Zheng Weihong, V.N. Kotov, and J. Oitmaa, 
unpublished. 

R.R.P. Singh and D.A. Huse, Phys. Rev. Lett., 68, 1766 
(1992). 

P. Lecheminant et al, Phys. Rev. B56, 2521 (1997). 
N. Elstner and A. P. Young, Phys. Rev. B50, 6871 (1994). 
R. R. P. Singh, O. A. Starykh and P. J. Freitas, J. Appl. 
Phys. 83, 7387 (1998). 

P. Chandra, P. Coleman, and A. I. Larkin, J. Phys.: Con- 
dens. Matter 2, 7 933 (1990). 

R. H. McKenzie, cond-mat/9802198, to appear in Com- 
ments on Condensed Matter Physics. 
For an example of uniaxial experiments on these materi- 
als see, C. E. Campos et al, Physica B 211, 293 (1995). 
J. Schmalian, Ph ys. Rev. Lett. 81, 4 232 (1998); H. Kino 
and H. Kontani, |cond-mat/9807147t H. Kondo and T. 



[28] H.A. Bethe, Z. Phys. 71, 205 (1931). 



Moriya, sond-mat/9807322; M. Votja and E. Dagotto, 



cond-mat /9807168|, to appear in Phys. Rev. B; K. Kuroki 



and H. Aoki, cond-mat/9812026 



S. R. White and I. Affleck, Phys. Rev. B 54, 9862 (1996), 
and references therein. 

R. R. P. Singh and N. Elstner, Phys. Rev. Let t. 81, 



4732 (1998); R. R. P Singh and Z. Weihong, |cond 



mat/9811028 



[15] A. V. Chubukov, S. Sachdev, and T. Senthil, Nucl. Phys. 

B 426, 601 (1994). 
[16] A. Angelucci, Phys. Rev. B 44, 6849 (1991); ibid. 45, 

5387 (1992). 

[17] N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 
(1991). 

[18] H.X. He, C.J. Hamer and J. Oitmaa, J. Phys. A 23, 1775 
(1990). 

[19] M. P. Gelfand, R.R.P. Singh, and D.A. Huse, J. Stat. 
Phys. 59, 1093 (1990). 

[20] M. P. Gelfand, Solid State Commun. 98, 11 (1996). 

[21] D.A. Huse, Phys. Rev. B 37, 2380 (1988). 

[22] W.H. Zheng, J. Oitmaa and C.J. Hamer, Phys. Rev. B 
43, 8321 (1991). 

[23] A.J. Guttmann, in "Phase Transitions and Critical Phe- 
nomena", Vol. 13 ed. C. Domb and J. Lebowitz (New 
York, Academic, 1989). 

[24] J. Merino, unpublished. 

[25] H. Kontani, M. E. Zhitomirsky, and K. Ueda, J. Phys. 

Soc. Jap. 65, 1566 (1996). 
[26] J. C. Le Guillou and J. Zinn- Justin, Phys. Rev. Lett. 39, 

95 (1977). 

[27] P. Azaria, B. Delamotte, and T. Jolicoeur, Phys. Rev. 
Lett. 64, 3175 (1990); P. Azaria, B. Delamotte, and D. 
Mouhanna, ibid. 68, 1762 (1992). 



11 



